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The consequences of spin-rotation-gravity coupling are worked out for linear gravitational waves. 
The coupling of helicity of the wave with the rotation of a gravitational- wave antenna is investigated 
and the resulting modifications in the Doppler effect and aberration are pointed out for incident 
high-frequency gravitational radiation. Extending these results to the case of a gravitomagnetic field 
via the gravitational Larmor theorem, the rotation of linear polarization of gravitational radiation 
propagating in the field of a rotating mass is studied. It is shown that in this case the linear 
polarization state rotates by twice the Skrotskii angle as a consequence of the spin-2 character of 
linear gravitational waves. 
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I. INTRODUCTION 



In a recent paper on the purely gravitational spin- 
rotation coupling, Shen Q has treated the coupling of 
graviton spin to the gravitomagnetic field. In this way, 
spin-gravity coupling has been extended to include gravi- 
tational waves. The subject of spin-rotation-gravity cou- 
pling for a particle of spin s has been reviewed in Refs. 
12 and and discussions of more recent advances are 
contained in Refs. [4-12]; however, these treatments have 
ignored the s — 2 case. Shen's field-theoretical approach 
is based on a weak-field approximation scheme that em- 
phasizes the self-interaction of the nonlinear gravitational 
field [2. 

The purpose of the present paper is to investigate the 
consequences of the helicity-gravitomagnetic field cou- 
pling for weak gravitational waves. To provide a compre- 
hensive treatment, we begin with the analysis of the prop- 
agation of free gravitational waves in a Minkowski space- 
time background from the standpoint of a uniformly ro- 
tating observer. We then generalize our helicity-rotation 
coupling results to the propagation of gravitational waves 
in the field of a rotating astronomical mass via the grav- 
itational Larmor theorem. 

The plan of this paper is as follows. In Section II, we 
study the reception of a free gravitational wave by a ro- 
tating observer in Minkowski spacetime. Section III deals 
with the modification of Doppler effect and aberration for 
gravitational waves caused by the helicity-rotation cou- 
pling. The results are extended to the gravitational case 
in Section IV using the gravitational Larmor theorem. 
It is then possible to study the influence of the gravit- 
omagnetic field of a rotating source on the propagation 
of high-frequency gravitational radiation. The rotation 
of the linear polarization state of gravitational radiation 
propagating in the field of a rotating source is directly 
calculated in Section V using the eikonal approximation. 
Section VI contains a discussion of our results. In this 
paper, we choose units such that c — 1, moreover, the 



signature of the metric is +2 in our convention. 



II. HELICITY-ROTATION COUPLING 

Imagine a class of uniformly rotating observers O' in 
a global inertial frame of reference. For the sake of con- 
creteness, we choose Cartesian coordinates such that the 
observers rotate with a frequency 17 about the z axis, 
each on a circle of radius p, 0<p<l/fi, parallel to the 
(x, y) plane of the inertial system. The local orthonormal 
tetrad frame of each observer is given in the {t,x,y,z) 
system by 



AM 
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7(1, —w sin 0,1) COS 0,0), (2.1) 

(0, cos 0, sin 0,0), (2.2) 

7(w, — sin0, COS0, 0), (2.3) 

(0,0,0,1), (2.4) 



where cj) — fit — jflr, v — Qp, 7 is the Lorentz factor 
corresponding to v and r is the proper time such that t = 
at t = 0. Regarding the motion of each observer, we 
note that Aj^^^^ and A'^^s^ indicate the radial and tangential 
directions, respectively, in cylindrical coordinates. 

Let us first consider an incident plane monochromatic 
gravitational wave of frequency lo and wave vector k = 
w(0,0, 1), so that each observer rotates about the direc- 
tion of wave propagation. The gravitational potential of 
the incident radiation is given by the symmetric tensor 
/ipjy, which represents a small perturbation of the back- 
ground Minkowski metric 77^^. Therefore, only terms 
that are linear in h^i, will be considered throughout. In 
the transverse-traceless gauge, /iqm = and the potential 
for circularly polarized gravitational radiation is given by 
the matrix (hij) = Re(P±), where 



P± = (ee ± «e®)^(k) 



(2.5) 



The upper (lower) sign corresponds to positive (nega- 
tive) helicity and the two independent linear polarization 
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states are denoted by 





(2.6) 



All of the operations involving wave functions arc linear; 
therefore, only the real part of the relevant quantities will 
be of any physical significance. 

Consider the measurement of spacetime curvature by 
the static observers O; the Riemann tensor can be ex- 
pressed as a 6 X 6 matrix (TZab), where the indices A and 
B range over the set {01, 02, 03, 23, 31, 12}. The gravita- 
tional field is given by the measured components of the 
Riemann tensor and for a Ricci-flat field. 



n = 



E H 

H -E 



(2.7) 



where E and H are 3x3 symmetric and traceless ma- 
trices corresponding to the electric and magnetic com- 
ponents of the curvature tensor. It turns out that for 
all of the circularly polarized gravitational waves consid- 
ered in this section, the gauge-invariant Riemann tensor 
is of the form of Eq. H2.7|l with H = ±iE, which, tak- 
ing due account of the proper definitions of curvature- 
based gravitoelectric and gravitomagnetic fields Ts'l , cor- 
responds exactly to the analogous electromagnetic case. 
It is therefore sufficient to focus attention only on the 
electric part of the Riemann tensor. For the static ob- 
servers O, the field of circularly polarized radiation is 
thus given hy E — ^uj^P±. 

The gravitational field as determined by the rotating 
observers C is given by i?^j,p^A^^^A||'^^A|'^^A^_5j. These 

may be expressed via the real part of TZ' = CIZC^ , where 
£ is a real 6x6 matrix that can be determined from Eqs. 
(|231)-(E!ill- We find that 



C = 



A B 

-B A, 



(2.8) 



where A and B are given by 







7COS0 7sin0 \ / ^ 

A = I -sin0 COS0 ,6 = W7 
7 / \ COS0 sini/) 



It is then possible to express TV as 

7^' = 



C± ±iC± 
±iC± -C± 



where C± is given by 
C± = -uj"^ ±ij -1 -W7 I h(k). 



9 9 9 

vj — w 7 



(2.11) 



Thus the frequency measured by the rotating observers 
via the temporal dependence of Eq. (|2.11() is 



(2.12) 



which clearly exhibits the contribution of helicity- 
rotation coupling {14] . This expression is the exact spin-2 
analog of the electromagnetic result |^, that has been 
observationally verified for a; ^ fi. Let us note that 
a simple application of the Doppler formula would lead 
to the transverse Doppler frequency oj'j^ = juj; however, 
the Doppler formula must be modified by taking into ac- 
count the helicity-rotation coupling as in Eq. (|2.12|) . For 
a packet of free gravitational waves propagating along 
the axis of rotation of the observer, the frequency of each 
Fourier component would be affected as in Eq. H2.12|l . 
A complete discussion of the physical implications of Eq. 
(|2.12|) will not be given here, since such treatment would 
be entirely analogous to the electromagnetic case that 
has been discussed in detail in Ref . . 

Let us now consider the extension of helicity-rotation 
coupling to the case of oblique incidence. Expressing the 
incident plane wave in terms of spherical waves whose 
dependence upon time t and the azimuthal coordinate 
(p is of the form exp(— iwf -I- irrup), and taking into ac- 
count the fact that a transformation to the frame of the 
rotating observer involves (r, i9, Lp) (r, Lp') in terms 
of spherical polar coordinates such that ip = Lp' + Qt, we 
find that 



j' = j{uj-mn), TO = 0,±1,±2,. 



(2.13) 



Here m is the multipole parameter such that mh is the 
total (orbital plus spin) angular momentum of the ra- 
diation field along the direction of rotation of the ob- 
server. When this direction coincides with the direction 
of wave propagation, only the spin contributes to m in 
Eq. (|2.13(l and hence we recover Eq. H2.12|l for radia- 
tion of definite helicity (m = ±2). It proves interesting 
to consider the other special case where the contribu- 
tion of the orbital angular momentum of the obliquely- 
incident radiation field vanishes, namely, the reception 
of the radiation by a rotating observer Oq at x = 0. In 
this case, the noninertial observer is at rest at the ori- 
gin of spatial coordinates, but refers its measurements 
to axes rotating with frequency fl. The tetrad frame 
for O'q is given by Eqs. H2.1|l - H2.4|l with v — and 
7 = 1. Imagine therefore an incident plane circularly- 
polarized gravitational wave with wave vector k = wk, 
where k = (0, — sin0, cos ^). According to the static in- 
ertial observers, the natural orthonormal triad for the 
wave is (x, fi, k), where fi = k x x = (0, cos 6*, sin 6*). It 
is straightforward to express the potential of the wave in 
the transverse-traceless gauge as (/ly ) = Re(P±), where 



P± = (e® ± «e®)/i(k)e' 



z(k . x—ut) 



(2.14) 



and 





' sin 9 cos ( 
- sin^ 6 
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cos 6 sin 9 





(2.15) 



The gravitational field measured by the static inertial 
observers, TZ, has the standard form and its electric part 
is given by \ui'^P±. Similarly, the field according to the 
special rotating observer Oq is TV with electric compo- 
nents given by ^oj^P'^. To determine Pj^, let 



cos 4> sin (p 
=1 — sin (/) cos (j) 
1 



(2.16) 



be the rotation matrix that relates the spatial frame 
(x',y',z') of the rotating observer to that of the static 
inertial observers. We find that 



P± 



RnP±Rl (2.17) 
It then follows after some algebra that is given by 



i{uj — 7nil)t 



(2.18) 
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±i 
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±i 


-1 
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±i 


-1 
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so that in this case the measured frequencies are Wg = 
CO - mn, where m = 0,±1,±2. Here TW^™) are 3x3 
symmetric and traceless matrices given by Al*^"^ — 
diag(i, ^,-1), 



moreover, the coefficients fi„i can be expressed as 
fiQ — sm^6, iii—±a±sm9, — {a-±Y 

where 

1 



a± = -(1 ± cos( 



(2.19) 

(2.20) 
(2.21) 

(2.22) 



It follows from these results that Oq can express its mea- 
surements as a Fourier sum of frequencies uj'q — lo — mfl 
with amplitudes given in Eqs. (|2.18|I - H2.22|I . It is rea- 
sonable to assume that the weight Wm, or intensity, as- 
signed to each Fourier component is given by the sum of 
the squares of the absolute magnitudes of the elements 
of the corresponding matrix. Thus, W,n — \h{'k)\'^Wm, 
where 



wq — - sm 
2 



Wi 
2 „• 2 n 



4(a±)2 sin^ e, W2 = 4(a±)^ (2.23) 
w-i = 4(a^)^ sin^ 61, w_2 = 4(azp)^ (2.24) 
and one can show that 

J2wm=4.. (2.25) 



Let us now define the relative weight of each frequency 
{lu — jnil) in the Fourier sum to be 



W„ 



1 



(2.26) 



therefore, the average frequency measured by the ob- 
server O'q can be computed and the result is 

(luq) ='^{lu - mfl)p,n = uj cos 9. (2.27) 



These considerations can be given a proper physical 
interpretation based upon the eigenstates of a particle 
with spin 2h. That is, according to the representation 
theory of the rotation group, the eigenstates of the par- 
ticle with respect to the coordinate system (x, fi, k) can 
be transformed to the (x, y,z) system using the matrix 

(■^mm') ^'^^ i = 2 Each element of this matrix is 

2. This 



given, up to a phase factor, by d^J} ,(9) with j 
latter matrix can be expressed [13 as 



/al 
bo 



-b+ 

c+ 
-Co 

c_ 



bo 

Co 

ao 
-Co 
bo 



c_ 

Co 

c+ 
b+ 



at \ 
-6_ * 

bo 

-b+ 

-.2 



(2.28) 



a. 



+ 



where a+ and a_ are given by Eq. 12.22|) and 

1 \/6 
flo = -(l+3cos26'), 60 = sin^ 0, h± — —a±sa\t 



Co = ^ sin 29, 



c± 



(2.29) 



- (cos 61 ± cos 26*). (2.30) 



Consider an incident graviton with definite helicity ±2?i 
as in Eq. (|2.14|) . The probability amplitude that the 
graviton has an angular momentum mTi along the z axis, 
corresponding to a frequency oJq — oj — mf2 as measured 
by O'q, is given up to a phase factor by Cm, where 



(Cm) 



/ 4 \ 

bo 

V at J 



(2.31) 



Here the upper (lower) sign refers to an initial incident 
state of positive (negative) helicity. Equation (|2.31|l is 
obtained from the first and last columns of the matrix 
(|2.28l) . since the state of the particle in the (x,y,z) sys- 
tem is obtained from the application of the matrix (|2.28() 
on the helicity states of the incident graviton. It follows 
that the probability that an incident graviton of helicity 
±2?i has spin mfi, m = 0,±1,±2, along the direction of 
rotation of the observer is given by |CmP, where 



Pr, 



(2.32) 
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based upon the comparison of Eqs. H2.26|l and H2.31|l . 
Thus pm is indeed the probabihty that observer Oq would 
measure frequency uj — mfi and hence the average mea- 
sured frequency is in fact (wg) given by Eq. (|2.27|) . 

Introducing the hehcity vector H — ±k, the average 
frequency measured by the observer can be written as 
(wg) = uj — 2fl . n, which is another expression of hehcity- 
rotation couphng. In view of Eq. (|2.12() . we may in- 
terpret the expression for (cjg) as follows: the rotation 
frequency of the observer SI = Qz, may be decomposed 
into a component of magnitude Q cos 9 parallel to the 
wave vector k and a component of magnitude sin 9 per- 
pendicular to k. On the average, the latter component 
does not contribute to the measured frequency; hence 
(wg) = u! ^ 2(ricos0) in agreement with Eq. (|2.12|l . 

The average frequency measured by the observer is 
expected to be the same as the result that would be 
obtained in the JWKB regime in accordance with the 
quasi-classical approximation. In the case under consid- 
eration, this corresponds to the high-frequency regime for 
gravitational waves; that is, waves with uj fl. To see 
how this comes about explicitly, we follow the approach 
that has been developed for the electromagnetic case Q 
and adapt it to the gravitational case under considera- 
tion here. It is first necessary to define a triad (a,/3, k) 
that can provide a more natural polarization basis for the 
rotating observer Og such that a and $ remain "fixed" 
in the rotating frame as much as possible. These unit 
vectors are defined by 



a — xcos$ + nsin$, 
P = — xsin (& -|- ncos 3>, 



(2.33) 
(2.34) 



where $ is given by 



sin<I> = — cos6'sinilt, cos<I> = — cosilt (2.35) 
D D ^ ' 



and D > can be obtained from 

= cos^ 9 + sin^ 9 cos^ ?Lt. 



(2.36) 



We note that $ reduces to fit for 6* = and to — ili for 
6* = TT, while $ = for 6* = 7r/2. 

The orthonormal triad (a, /3, k) is related to (x, y, z) by 
a rotation; in fact, the transformation that takes (x, y, z) 
to (a, /3, k) is given by 



cos $ cos 9 sin $ sin 9 sin $ 
T = ( — sin $ cos 9 cos sin 9 cos $ 
— sin 9 cos 9 



(2.37) 



We wish to find the new polarization basis for linearly po- 
larized gravitational waves based on the new triad; this 
can be achieved by a similarity transformation of the ba- 
sis given in Eq. (|2.6|l . and the result is 



It is possible to connect this new basis with Eq. (|2.15f) 
as follows 

i = ?eCOs2$-|-?55sin2$, (2.39) 

V = -eg, sin 2$ -I- COS 2$. (2.40) 

The rotation by 2$ takes a simple form for the circular 
polarization basis, namely, 

e±iz^= (e~©±^e«)eT2.*^ (2.41) 
Using this relation in Eq. H2.14|l . we have 



P± = (e±ij^)/i(k)e-*"*±"**, 
so that from Eq. (|TT7|) . 



Pi = (e' ± w')hiy.)e 



where the new polarization basis involving 



(2.42) 



(2.43) 



(2.44) 



is naturally adapted to the rotating observer in the sense 
that the temporal dependence in Eq. (|2.43() has been 
transferred to the phase of the wave as much as possible. 
Let us note that ^' and v' are obtained from and e^, 
respectively, by a unitary transformation involving the 
orthogonal matrix RqT'^ . 

According to Eq. H2.43|l . observer Oq receives a cir- 
cularly polarized wave with a phase —cut ± 2$. The fre- 
quency of the wave is defined to be the negative gra- 
dient of the phase with respect to time; hence, Wg — 
UJ T 2d<^/dt. It follows from Eq. ifT^ that 



1)1 



VL cos 9 



(2.45) 



In practice, the frequency determination would necessi- 
tate the reception of at least a few oscillations of the 
wave. During such a period of time t, assuming that the 
observations begin at t = 0, e' = fit ^ 1 due to the fact 
that ri ^ w; moreover. 



(2.46) 



(2.38) 



by Eq. Thus uj'q = uj ^ 2^008 9 in the high- 

frequency regime. This result is equivalent to the intu- 
itive expectation that from the standpoint of the rotating 
observer, the spin of the radiation field should precess in 
the opposite sense. This circumstance can be restated in 
terms of the rotation of the state of linear polarization of 
the gravitational wave as explained in Section IV. 

Henceforward, we limit our considerations to the high- 
frequency regime {oj ^ il). It is then possible to gener- 
alize the main result that we have obtained for the fixed 
rotating observers Oq to the case of arbitrary rotating 
observers and thereby determine the modifications in the 
Doppler effect and aberration that are brought about by 
the hclicity-rotation coupling. This is done in the next 
section. 
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III. MODIFIED DOPPLER AND ABERRATION 
FORMULAS FOR GRAVITATIONAL WAVES 

In a background global inertial frame, we first consider 
the class of noninertial observers that are at rest in the 
inertial frame but refer their measurements to rotating 
axes as discussed in Section II. Each of these observers 
carries a tetrad frame of the form 

A'(^) = (1,0,0,0), A;^) = (O,cos0,sin0,O), (3.1) 

A'(^) = (O,-sin0,cos^,O), A;^) = (0,0,0,1), (3.2) 

where (j) = ilt a.s before. We are interested in the deter- 
mination of frequency and wave vector of incident high- 
frequency (w ^ f2) gravitational waves. To this end, 
we focus attention on the noninertial observer O'q at the 
origin of spatial coordinates and imagine the class of ob- 
servers that are at the rest in the rotating system of Oq, 
each with a tetrad of the form of Eqs. H2.1|l - H2.4|l . The 
gravitational field measured by this class of rotating ob- 
servers is given by TZ' = CR.C'^ . A complete Fourier 
analysis of TV in time and space is needed to determine 
the frequency and wave-vector content of the incident 
radiation according to Op . In the high-frequency regime 
{uj 3> fi), the measurements of Oq can be restricted in 
space to the cylindrical domain of radius <C Vl^^ . The 
analysis of frequency determination for lo ^ ^ has been 
given in Section II and a corresponding analysis of the 
wave-vector determination will not be carried out here, 
since it is entirely analogous to the electromagnetic case 
presented in detail in Ref. It follows from an analy- 
sis similar to the one given in Section 2 of Ref. (l^ that 
kg = ko. Thus we conclude that to lowest order 

ujQ = uj-sfi.n, k[, = k (3.3) 

for Og in the high-frequency regime. The dispersion re- 
lation for Og is then Wg = fcg =F skg . f2, where fcg = |kg|. 
Here s = 1 for electromagnetic waves and s — 2 for grav- 
itational waves. Indeed, these results hold for each mem- 
ber of the rotating class of observers that are at rest and 
carry the tetrad frame A|^j given by Eqs. 1)3. 111 - 13. 2|l . 

The generalization of Eq. 1)3. 3|l to the case of rotating 
observers that are not at rest in the background global 
inertial frame can be simply obtained from the obser- 
vation that at each event along the circular path of an 
observer O' , its tetrad A^^ > is related to the tetrad A'.''^ 
of the static rotating observer at that event by a Lorentz 
boost 

Aro)-7[A;i;)+«A;^)], Af^^^A;';), (3.4) 

Af,)^7[A;^) + -A;^;)], a|'3)^a;^;). (3.5) 

It follows from Lorentz invariance that (w', k') for O' are 
related to (a;g,kQ) by the standard Doppler and aberra- 
tion formulas; that is 

Lo' = 7(4-v.k;,), (3.6) 
k' = k[, + ^l^(v.k;,)v-74v. (3.7) 



Substituting Eq. in Eqs. (|^^ - (|T7|l . we obtain the 

modified Doppler and aberration formulas for gravita- 
tional waves in the high-frequency regime 

uj' = 7[(cj-sH.O)-v.k], (3.8) 
k' = k+-^^^^^(v.k)v-7(w-sH.n)v. (3.9) 

V 

Eq. H3.8(l may be interpreted in terms of Eq. (|2.13() in 
the eikonal approximation, namely, 

tj' = 7(cj-j.rj), j = rxk-hsH, (3.10) 

where ?ij is the total angular momentum of the graviton 
(s = 2) or the photon {s = 1). The results of this section 
for 5 = 2 are expected to be of importance in experi- 
ments involving the reception of gravitational waves by 
antennas rotating with frequency f2 <C w. Clearly, cur- 
rent large-scale Earth-fixed gravitational-wave antennas 
rotate with the frequency of rotation of the Earth. 



IV. HELICITY- GRAVITY COUPLING 

It is possible to employ Einstein's principle of equiv- 
alence in order to extend the helicity-rotation coupling 
discussed in previous sections to the propagation of grav- 
itational waves in a curved spacetime background. For 
this purpose, we consider a solution of the linearized 
gravitational field equations. We assume that this back- 
ground field is due to isolated gravitating sources that 
move slowly compared to the speed of light in vacuum. It 
is possible to describe such a background field in terms of 
gravitoelectromagnetism ( "GEM" ) in close analogy with 
Maxwell's electrodynamics. GEM has been thoroughly 
reviewed in Ref. and references cited therein. In this 
approach to GEM, the metric of the curved spacetime 
background is of the form 

- (1 - 2f)dt^ - 4(S . dx)rft + (1 + 2f)5,jdx^dx^ , (4.1) 

where in our convention / and S are respectively the 
gravitoelectric and gravitomagnetic potentials subject to 
a ("Lorentz") gauge condition 

| + V.(is)=0. (4.2) 

These potentials are connected to sources via the lin- 
earized gravitational field equations. For the treatment 
in this section and in conformity with our general linear 
approach, we assume that in the presence of an incident 
gravitational wave the deviation from the Minkowski 
spacetime is a linear superposition of the perturbations 
due to the isolated sources and the incident wave. 
The GEM fields are defined by 

1 

F = -V/-- — , B-VxS, (4.3) 
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in our special convention that has been designed to 
provide the closest possible connection with the stan- 
dard formulas of classical electrodynamics. Eqs. (|4.2|) 
and H4.3|l together with the linearized gravitational field 
equations lead to the Maxwell equations for the GEM 
fields O. 

It follows from a detailed discussion of Einstein's prin- 
ciple of equivalence within the context of GEM that the 
gravitoelectric field is in effect locally equivalent to a 
translationally accelerated system, while the gravitomag- 
netic field is in effect locally equivalent to a rotating sys- 
tem Traditionally, Einstein's heuristic principle of 
equivalence refers to the accelerated "elevator" in rela- 
tion with the gravitoelectric field of the source; however, 
the rotation of the elevator is in general necessary as well 
to take due account of the corresponding gravitomagnetic 
field. 

In keeping with the electromagnetic analogy, this ap- 
plication of Einstein's principle of equivalence is the con- 
tent of the gravitational Larmor theorem [l^ IT^ . It 
turns out that a spinning particle at rest in the exte- 
rior field of a rotating mass precesses with frequency B, 
which is equivalent to what would be observed from a lo- 
cal frame of reference rotating with frequency $7 = — B. 
Following this general line of thought, we may conclude 
from the results of Section III that for observers at rest 
in the exterior GEM background, the local dispersion re- 
lation 



UJ 



/cisk.B 



(4.4) 



is approximately valid for high-frequency incident grav- 
itational (s — 2) or electromagnetic (s = 1) waves. Eq. 
(|4.4|l follows from Eq. H3.3(l with Q, -B in accor- 
dance with the gravitational Larmor theorem. Recog- 
nizing that the local dispersion relation H4.4|l ignores the 
usual global GEM effects, such as the bending of the in- 
cident beam of radiation, one may nevertheless employ 
Eq. H4.4|l globally in order to uncover the specific con- 
sequences of helicity-gravity coupling. The results may 
then be superposed on the standard GEM effects in line 
with our general linear perturbative approach. 

An interesting consequence of the coupling of helicity 
with the gravitomagnetic field is the rotation of the state 
of linear polarization of a gravitational wave that prop- 
agates in the field of a rotating mass. To illustrate this 
effect, we assume that the background field is stationary 
and is due to a rotating astronomical source. The GEM 
fields have been determined in this case [3|-The curl of 
the gravitomagnetic field B vanishes in the exterior of 
the source; hence. 



B 



-VQ, 



(4.5) 



where Q is the gravitomagnetic scalar potential. Far from 
the source 



/ 



GM 
r 

GJ.r 



GJ X r 



C T 

B^^[3(J.f)f-J], 



(4.6) 
(4.7) 



where M is the mass and J = Jz is the angular momen- 
tum of the source. 

Consider a linearly polarized gravitational wave start- 
ing at z — zq far from the source and propagating out- 
ward along its rotation axis. Let II^"' and IIj"' be the 
linear polarization tensors for the wave. Assuming that 
at z = Zq the state of the wave is given by the real part of 
tp^^'^ = ipH^'^ exp{—iujt), where ip, <C 1, is a constant 
amplitude, then for any z 



Here S is given by 



S± = k±dz, 



(4.8) 



(4.9) 



where k+ and are the wave numbers of the posi- 
tive and negative-helicity components of the gravitational 
wave. Expressing S± as 

S+^Sq~ a, S-^So + a, (4.10) 

we find that 4'^" can be written as 

V'^'^ V;(n5"'cosA-hn^''sinA)e*'^°-*"*. (4.11) 

Inspection of this equation reveals that as the wave prop- 
agates, the linear polarization state rotates by an angle 
A given by 



A.i 
2 



(fc_ — k-^)dz. 



(4.12) 



To compute this angle, we write Eq. (|4.4|l for the positive 
and negative helicity components of the wave, 

UJ = k+ + sB^, Lu^k^-sB,^, (4.13) 



so that we find 
A = s 



B,dz ^ s[Q{zo) - Q{z)]. 



Using Eq. H4.7|) . we finally have 



A = sGJ{ 



1 



(4.14) 



(4.15) 



For electromagnetic radiation (s = 1), this gravitomag- 
netic rotation of the plane of polarization was first stud- 
ied by Skrotskii 0|; detailed discussions and references 
are contained in Ref. |23| and references cited therein. 
The angle of rotation for gravitational radiation (s — 2) 
is twice the Skrotskii angle. Let us note that A van- 
ishes for waves propagating from — oo to -|-oo along the 
z axis; this is an instance of a general result discussed 
in Appendix A, where the consequences of Eq. 14.4|l are 
worked out in a more general context. 

The angle of rotation of the state of linear polarization 
A is independent of the frequency (or wavelength) of the 
radiation; therefore, Eq. H4.15(l is valid in the limit of 
vanishing wavelength, namely, the JWKB (or eikonal) 
limit. In this limit, gravitational waves propagate along 
a null geodesic; in this full treatment is contained 

in the next section. 
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V. ROTATION OF LINEAR POLARIZATION 

The purpose of this section is to compute (within the 
eikonal approximation scheme) the rate of rotation of 
the state of hnear polarization of high-frequency grav- 
itational waves propagating in the exterior background 
field of a rotating astronomical source. 

Consider the propagation of gravitational radiation on 
a background spacetime such that the waves cause a 
small perturbation. Let g^i^{x) be the metric tensor of 
the background field in a given coordinate system and 
giiu = giiu-^h^u be the spacetime metric tensor. Under an 
infinitesimal coordinate transformation x'^ = x'^ — e^{x), 

h'^Ax) ^ h^^{x) + e^\^, (5.1) 

where the vertical bar denotes covariant differentation 
with respect to g^„. We also raise and lower indices, etc., 
with gp^u. Thus the perturbation is determined up to a 
gauge transformation given by Eq. H5.1|l . Introducing 
the trace-reversed potential 

i't^u = h^,^ - ^g^^g^^hpa, (5.2) 
we find that under a gauge transformation 

"^^'t^u = ''Pi^u + e^ii. + Ci/Ip - g^,u (5.3) 

It is convenient to impose the transverse gauge condition 

V-^V = 0, (5.4) 

which does not fix the gauge completely, however. It 
turns out that any solution of 

e^l^" -t- W^e" = (5.5) 

in Eq. leads to V'''^"],, = if Eq. is assumed. 

In the transverse gauge, the gravitational field equations 
imply that V'/^i/ satisfies the wave equation plj 

V-^,,/ + 2R^p,,^P" = 0, (5.6) 

where is assumed to be Ricci-fiat in the spacetime 
region under consideration here. 

To describe the propagation of the wave function 
in the eikonal approximation, we seek a solution of Eq. 
(|5.6|) in the form 

V',.. = Re{x^,(x;e)e'^^"''^(^)}, (5.7) 

where e, < e <C 1, is directly proportional to the wave- 
length of the radiation. In the eikonal approximation, 
X^i,(a;;e) is expressed as an asymptotic series in powers 
of e 

Xtiv^^x; e) = Xy^v{x) + tp^,y{x) + (^K^,y{x) + .... (5.8) 

Let fc^ = da{x)/dx'^'' be the propagation vector of the 
wave; then, the substitution of equations l|5.7|) and (|5.8|) 



in the gauge condition l|5.4|) and propagation equation 
(|5.6() results in series that contain powers of 1/e in addi- 
tion to powers of e. It follows from Eq. (|5.4|) that there 
is only one such term involving 1/e and in the eikonal 
limit (e — > 0), the coefficient of this term must vanish; 
therefore, 

X^.uk'' = 0. (5.9) 

Moreover, the propagation equation 1)5. 6|l involves 1/e^ 
and 1/e terms and the coefficients of these terms must 
also vanish in the eikonal limit, hence we have respec- 
tively 

k^k^" = 0, (5.10) 

and 

*(2x,.H-fc" + Xm-^V) - = 0. (5.11) 

Let us first note that equation H5.10|l implies that the 
radiation follows a null geodesic in the eikonal limit. It 
follows from fc^, = da{x)/dx^ that fc^|^ = fc^i^. Thus, 
taking covariant derivative of Eq. H5.1U|I . we get that 
kfj.\uk'^ = 0; hence, 

= 0. (5.12) 

The geodesic equation follows from Eq. H5.12|l and fc^ = 
dx'^/dX, where A is an affine parameter along the path. 

Equations H5.9|l and (|5.11|l describe the propagation of 
the wave amplitude along the null geodesic, since it 
follows from Eq. H5.11|l that 

2x^..\ak^ + Xt..k^l, = 0. (5.13) 

An immediate consequence of this relation is that S*^ = 
X^^X^"^ satisfies the conservation law 

(S°fc-)|, =0, (5.14) 

which can be interpreted as the conservation of the 
"graviton" number along the null geodesic congruence. 
These results have been based on terms involving 1/e 
and 1/e^; taking account of the other terms in the eikonal 
series, i.e., those involving e",n = 0,1,2,..., would sim- 
ply specify the manner in which the general wave ampli- 
tude Xfj_iy(x; e) varies along the null geodesic. This eikonal 
(or JWKB) treatment of gravitational radiation has been 
previously considered in Ref. Q^. A critical assessment 
of the eikonal approximation scheme is contained in Ref. 

In the eikonal approximation scheme, the curves x'^ — 
x^{X) that have k^ — dx^ /dX as tangent vectors are null 
geodesies orthogonal to the surfaces of constant phase cr. 
Imagine a bundle of such null rays in a congruence char- 
acterized by the propagation vector k^. Let us define a 
null tetrad system {k^ ,1^ ,n*^) such that k'^lf_t = -1 
and n^n*^ = 1 are the only nonvanishing scalar products 
among the four null vectors. Starting from an observer's 
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orthonormal tetrad frame X'f x, the null frame is con- 

(a)' 

structed as follows: 
/c^ = -T^iKo) + Ks)^^ — Tk^Kq) ^ Ks)^ (5.15) 



V2' 



iV2 



are real, while n'^ and its complex conjugate n*'^ are com- 
plex, since is defined by 



~ 71 +'^('2)] 



(5.16) 



Here a = — v^fc^A^p^ is a nonzero constant. The tetrad 
system is assumed to be parallel propagated along the 
congruence. 

The symmetric and transverse tensor can be lo- 
cally expressed in terms of the parallel-propagated null 
tetrad as 

X^^ = ^+n^n^ + $-n*n* -|- k^T^ + k^T^. (5.17) 

It turns out that in this expansion n^rii/ corresponds to a 
positive helicity wave and n* corresponds to a negative 
helicity wave; thus, Xfj-v consists in general of a positive 
helicity part with amplitude , a negative helicity part 
with amplitude and a gauge part involving a trans- 
verse vector such fcT^ ~ (see the appendix of Ref. 
0). 

We note that (i/-^"^^^) = T.°/2 and 



so that the intensity of the wave depends only on its 
irreducible part 



(5.19) 



■0/^1/ = $+n^nj, -I- ^^n*nl. 



Using this irreducible part of x^^.v^ from Eq. (|5.13| 

that 



= 0, 



(5.20) 



where = (1/2)A:°'|^ is the expansion of the null con- 
gruence k'^. Let A be the area of the cross section of a 
bundle of null rays in the congruence; then, 



dA 



2M. 



(5.21) 



It follows that |<i>+pyl and |<i>_pA are separately con- 
served along the trajectory; that is, the number of pos- 
itive or negative helicity null rays ( "gravitons" ) is in- 
dependently conserved along the congruence. This ap- 
proach to the polarization of gravitational waves in the 
eikonal limit demonstrates that the parallel transport of 
polarization tensors n^ni^ and n* n* for positive and neg- 
ative helicities, respectively, is a natural interpretation of 
Eq. H5.13|) . This circumstance is a spin-2 analog of the 
spin-1 electromagnetic case, where the irreducible part 



of the wave amplitude corresponding to the vector po- 
tential is of the form (j)+n^ + (p-n*^. Furthermore, it is 
possible to define Stokes parameters E" for gravitational 
radiation such that r]a0T,°'T,^ — and 

E^ = E^ = -i($+$l-$_$;), (5.22) 



E-^ = 1$, K - |$_ 



(5.23) 



along the lines developed in Ref. |2J|. 

It is important to discuss the uniqueness of the repre- 
sentation H5.17|l . At any given event, the observer tetrad 
frame X'^^^ is unique up to a Lorentz transformation. We 
are interested in a subgroup of the Lorentz group that 
preserves k^; in fact, this is the little group of A:^ that is 
isomorphic to the Euclidean group in the plane. There- 
fore, under the action of the little group, fc'^ = k^, 



I'' = I'' + bn'' + b*n*'' + \b\^k^', 



-iS 



[n.^^ + b*k^'). 



(5.24) 
(5.25) 



Here 9 is real and corresponds to the constant angle of 
rotation in the {X^^iy X^^^-^) plane. For 9 = 0, Eqs. (|??^ 

and H5.25|l with a complex constant b correspond to the 
Abelian subgroup of the little group under which 



(5.26) 



(5.18) where /c^G^ — and G^ is given by 



Gp = $+6*71^ + + -($+6*2 + $_62)fc^. (5.27) 

Equation ()5.26(l amounts to a gauge transformation of 
^->^u\ therefore, under this gauge subgroup the Stokes pa- 
rameters for the radiation field remain invariant. On the 
other hand, with 6 = and under rotation of angle 9, 
^fj^u remains invariant if 



(5.28) 



which demonstrates that ($-) is the amplitude of the 
graviton in the positive (negative) helicity state. More- 
over, the Stokes parameters undergo a rotation as well, 
since a rotation of the angle 9 in the (A^^j,A^2)) plane 
induces a rotation of angle —49 in the (E^,E^) plane 

a. 

It proves useful to define real linear polarization tensors 
H^" and Hj" such that 



V2 



(Hi'^ + zH^"). 



(5.29) 



Here Hi and H2 are independent polarization states such 
that tp^^ can be written as 



(5.30) 
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where Li and L2 are the hnear polarization amphtudes. 
Under a constant rotation of angle in a plane per- 
pendicular to the spatial direction of propagation of the 
wave, 



n^"" cos 26 



av 



[26 



n^'' sin 26, 
n^'' cos 26, 



(5.31) 
(5.32) 



so that one linear polarization state turns into another 
under a rotation of 6 = 7r/4. In general, 6 — 7r/(2s) 
is the angle "between" the linear polarization states Hi 
and 112. 

The constant parameter a in the definition of fc'^ and 
in Eq. H5.15|l is related to the choice of the affine pa- 
rameter A, which is defined up to a linear transformation 
A ^ A' = constant + X/Aq where y4o ^ is a constant. 
Under this afhne transformation X'^^s^ is unchanged, but 

a aAo and hence fc'^ Aokf^, l^^ A^^lf^ and n'' is 
unchanged. A null rotation of the null tetrad is defined 
to be a 4-parameter group that consists of the combined 
action of the little group together with an affine trans- 
formation. It is the most general transformation that 
leaves the spatial direction of wave propagation vector 
k*^ invariant. 

It is important to note that under a rotation with an- 
gle 6, the linear polarization states of a massless spin-1 
field would rotate by 6, while that of a spin-2 field would 
rotate by 26. It has been shown in Ref. |23| that in a 
general gravitomagnetic field B, the plane of linear polar- 
ization of electromagnetic radiation rotates by an angle 



O-Skrotskii — / B . dx, 



(5.33) 



where the integral is evaluated along the spatial path of 
the null ray; see Section VII of Ref. |23| for a detailed 
derivation. The Skrotskii effect is the gravitomag- 
netic analog of the Faraday effect. Therefore, the state 
of linear polarization of gravitational radiation would ro- 
tate by an angle 2askrotskii in a general gravitomagnetic 
field. Thus, we recover, in the eikonal limit, the result of 
Section IV. 



VI. DISCUSSION 

It is useful to provide estimates of the spin-rotation- 
gravity coupling effects presented in this work. For 
Earth-based gravitational-wave antennas that rotate 
with the Earth, the effective rotation frequency is there- 
fore about 10~^Hz in Eqs. so that the in- 
cident gravitational waves that would be relevant in 
this case satisfy the high-frequency condition, namely, 

S> ri. Ignoring the helicity-rotation coupling would 
introduce a small systematic Doppler bias of magnitude 

2n/uj. 

Let us next consider, as in Section IV, a linearly polar- 
ized gravitational wave that propagates outward to in- 
finity starting from the north pole of a rotating astro- 
nomical system of radius Rq. We assume that uj ^ Bq, 



where Bq ~ 2GJ/Rq is the gravitomagnetic (Larmor) 
frequency of the system. The rotation of the linear po- 
larization state of the wave is in the same sense as the 
rotation of the source and the net angle of rotation is 
given by Eq. H4.15|l . namely, A = 2GJ/Rq. For a ho- 
mogeneous sphere of mass M rotating with frequency 17, 
A = O.SGAffi. This amounts to A « 0.025 rad (or about 
1.5°) for a millisecond pulsar of mass Ad ~ Mq and ro- 
tational period « lO^'^ s. For the Earth, however, the 
corresponding result would be negligibly small, that is, 
A « 10-15 rad. 



APPENDIX A 

The purpose of this appendix is to work out to first or- 
der in the helicity- gravity coupling the solution of Hamil- 
ton's equations for the dispersion relation (|4.4|) . The 
Hamiltonian for the ray motion is given in this case by 



n{r,k) = fc±sk.B(r) 



Hamilton's equations are 

dr 

di 
dk 

It 



-V[±sk.B(r)], 



where 



v, = |^ = i(k±.B) 
ok uj 



(Al) 

(A2) 
(A3) 

(A4) 



is the group velocity of the rays to first order in the cou- 
pling to the gravitomagnetic field. The background is 
stationary; therefore, for any ray that is a solution of 
equations l|A2|l and l|A3|l . uj = 7i(r,k) is a constant of 
the motion. 

Let (r+,k+) denote the solution of the equations of 
motion for a positive-helicity ray and (r_ , k_ ) denote the 
corresponding solution for a negative-helicity ray. Work- 
ing to first order in the helicity-gravity coupling, we let 



k+ = ko + K(<), k_ = ko — K(i), 
r+ = ro + q(i), r„ = Tq - q(t), 



where kg is constant, lo — |ko| and 
dro - 



(A5) 
(A6) 



(A7) 



is the equation of motion of the unperturbed ray in 
the absence of spin-gravity coupling. Substituting Eqs. 
(|X5l) - (|X7|) in the equations of motion (|X2|) - HA3p . we find 



= -{k + sB), 
dt u) 

^ = -.V(ko.B). 



(A8) 
(A9) 
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Using the fact that B = — VQ, we can write V(ko • B) = 
(ko . V)B, which can be expressed as dB/dt via Eq. (|A7|) 
in our approximation scheme. Hence it follows from Eq. 
(|A9|) that K+sB is a constant of the motion. This implies 
that the right-hand side of Eq. (|A8|) is a constant as well; 
therefore, 

K + sB = wVo, (AlO) 

where Vq = dq/ dt is a constant and B is evaluated along 
the average (unpolarized) ray. The requirement that lu 
be the same for both rays can be implemented using Eqs. 
(|A1I) and HA5|) and the result is 

ko .(k + sB) = 0, (All) 

so that Vg . kg = 1 . This completes the solution of the 
equations of motion. 

It follows from q(i) = q(0) + Vot and Eq. (jXeJ that 
the positive and negative helicity rays diverge away from 
the path of the average (unpolarized) radiation, since 
ko-Vo = 0. The integration constant Vq must be de- 
termined from the boundary conditions. Consider, for 
instance, the emission of rays of radiation that originate 
at z = zq on the axis of rotation of an astronomical mass. 
Let the initial direction of propagation of the radiation 
be normal to the z axis. Then, it follows from these ini- 
tial data that Vq = 2sG3 /{uJZg). Thus there will be a 
differential deflection of the radiation such that the posi- 
tive and negative helicity rays separate at a constant rate 
about the average direction of propagation kg . The total 
gravitomagnetic splitting angle between the rays would 
be 2Vo = 4:sGJ/{ujzq), which amounts to a few degrees 
for gravitational waves of frequency 10"^ Hz grazing the 
north pole of a millisecond pulsar. 

We emphasize that the polarization-dependent differ- 
ential deflection of rays will occur, in our treatment, only 
for radiation that originates near a rotating source. The 
gravitomagnetic field B rapidly falls off to zero far from 



the source; therefore, there will be no differential de- 
flection of rays in scattering situations. This conclusion 
agrees in the electromagnetic case (s = 1) with the re- 
sults of recent investigations 0|. 

One can use the results of this appendix to estimate 
the polarization-dependent time delay in the arrival of 
rays originating near a rotating astronomical source ^3 ■ 
Moreover, one can generalize the treatment of the rota- 
tion of linear polarization given in Section IV. In fact, it 
follows from Eqs. (|X5|) . HA10|) and (|A11|I that 

i(k_-k+).ko = sko.B, (A12) 

where ko . B — —dQ/dt along the unperturbed ray. 
Therefore, we find, following a method analogous to that 
used in Section IV, that the angle of rotation of linear 
polarization is 

A = s[Q{Yi) - Q{Yf)] (A13) 
for a linearly polarized ray of radiation traveling from 
Tj to ry. Thus for rays that are incident from infinity 
and travel to infinity, the net angle of rotation of linear 
polarization vanishes. 

Finally, we should mention that one can think of the 
effects discussed here in terms of the helicity dependence 
of the index of refraction of an effective medium for the 
propagation of rays. From the definition 

|k±|=c^n±, (AM) 

we find that 

n± = no T -k.B, (A15) 

where tiq = 1, since we have neglected here the 
polarization-independent bending of rays given by no ~ 
1 + 2/. 
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